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ABSTRACT 

We present the results of a high accuracy (cr ^ 0.005%) polarization monitor- 
ing of the Be Star Achernar that was carried out between July 7th and November 
5th, 2006. Our results indicate that, after a near quiescent phase from 1998 to 
2002, Achernar is presently in an active phase and has built a circumstellar 
disk. We detect variations both in the polarization level and position angle in 
timescales as short as one hour and as long as several weeks. Detailed modeling 
of the observed polarization strongly suggests that the short-term variations orig- 
inate from discrete mass ejection events which produce transient inhomogeneities 
in the inner disk. Long-term variations, on the other hand, can be explained by 
the formation of an inner ring following one or several mass ejection events. 



Subject headings: polarization 
{a Eridani) 



stars: emission line, Be — stars: individual 



1. Introduction 



A common feature to all Be stars is the presence of ionized circumstellar gas ( j Collins 



19871 ). Despite the general consensus that this gas is mainly distributed in a circumstellar 
disk, there is still much controversy on how the disk is built. It is generally belie ved that 
rotat ion plays an important role in the mass transfer between the star and the disk f Owocki 



but since evidence points to the fact that not all Be stars rotate subcriticallv fICranme: 



2005il. it follo ws that additional driving mechanisms must be present ( iTownsend et al 
Owockill2005l ). 
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This gap in our knowledge about disk formation can be attributed, at least in part, 
to the lack of sufficiently strong observational constraints. When the disk is fully formed, 
changes in its inner part may be obscured by the large emission from the outer parts. For 
this reason, the best phase to investigate the mass loss is when the disk has just started 
forming. 

Disks of Be stars show temporal variations on several timescales. Short- term (minutes 



to ho urs) line profile variations frequently occur in early-type Be stars (e.g. iPorri fc Stalio 



19881 ). The very short period of those variations suggests that they are associated with 
changes in the photospheric level or the innermost part of the disk and, thus, renders short- 
term variations as important properties to identify the additional mechanism(s) required for 
a rapidly rotating B star to become a Be star. 



Several such mechanisms were put forward (see lPorter fc Riviniusll2003l for an overview). 
Whatever the mechanism for ejecting gas from the photosphere is, it must involve the mass 
transfer from some part of the stellar photosphere to some position of the inner disk. It 
follows that if we can determine observationally the mass and geometry of the material as 
it is fed into the disk then we can begin to learn something about the physical mechanisms 
responsible for building the disk. 

In this letter we provide an attempt towards this goal. We present high-precision po- 
larimetric observations of the Be Star Achernar (HD10144, a Eri), with temporal resolutions 
of a few minutes. We achieved accuracies as high as 0.002% in the linear polarization and 
0.4° in position angle, depending on the night conditions. 



Observations 



The polarimetric observations have been performed with the 1.6-m Perkin-Elmer tele- 
scope and the 0.6-m Boiler & Chivens telescope at O PD/LNA - Brazil . We have used a 



CCD camera with a polarimetric module described in iMagalhaes et al. I (Il996l ). consisting 



of a rotating half-waveplate and a calcite prism placed in the telescope beam. A typical ob- 
servation consists of 8 consecutive wavepl ate positions (hereafter WPP) separated by 22?5. 
Details of data reduction can be found in IMagalhaes et al.l (119841 ). 



We define the average polarization of each WPP sequence, P, as the average of the M 
independent sets of 8 consecutive WPPs, i.e., P = M~^Ylj^Li Pj^ where Pi is the polarization 
for WPPs 1 to 8, Fa for WPPs 9 to 16, and so forth. Similarly, a = M'^ (E^|)^^^ 
summary of our polarimetric data is presented in Table [TJ 



A 
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In each observing run at least one polarized standard star was observed to calibrate the 
observed PA. The instrument stability was assessed by obs erving targets with a polarization 



level similar to Achernar. They were selected from the iHeiled (120001 ) catalog using the 
following criteria: (i) possess a polarization of about 0.15% known with good accuracy, and 
(ii) be of spectral type AV or later, to ensure that the polarization is purely interstellar. 
The selected targets were HD 150764, HD 118507 and HD 204867. We monitored those 
stars in the B filter for several hours, with a S/N better than 5000, thus producing several 
independent sets of 8 consecutive WPP positions. Our results indicate that the polarization 
is stable within a la level, even under less than ideal sky conditions. However, we found 
that the PA does vary more than la, but never more than 2a. From those tests we conclude 
that the instrument is stable, at a minimum, to the ~ 0.003% level for P and ~ 1.4° for the 
PA. 

Spectroscopic observations were carried out in November 05th, 2006 with the Coude 
spectrograph at the 1.60m telescope of OPD/LNA. We used a 600 groove/mm grating blazed 
at 6563 A in first order, resulting in a reciprocal dispersion of 0.24 A/pixel. The typical S/N 
ratio was 300. 



3. Results 

In Figure [1] we show the results of our polarization monitoring in the B filter between 
July 7th and November 5th, 2006. In one week, from July 7th to 14th, P increased from 
0.145 ± 0.003% to 0.188 ± 0.009%, a 30% increase. Two months later, by September 14th, 
the polarization had gone back to the 0.14% level and by November 5th it was as low as 
0.121 ± 0.002%. From the above, we can establish that the polarization level of Achernar 
has varied in the period of the observations with timescales of weeks to months. 

In addition to this long-term variation, the polarization seems to have varied also in 
much shorter timescales. To improve the temporal resolution of our measurements, we 
now group the WPPs in the following way: group 1, WPP from 1 to 8; group 2, WPP 
from 2 to 9; etc. This procedure effectively increases the temporal resolution, since two 
consecutive WPPs are separated of each other by about 4 minutes. However, this is not 
the true temporal resolution, because a polarimetric measurement is obtained from 8 WPPs 
and, thus, encompasses about l/2h of observations. 

Our results are shown in Figure [2|, where we plot the polarization vs. time for the 
WPP sequences of September 14th and November 5th. In September 14h, for instance, the 
polarization varied between 0.122 ±0.004% and 0.172 ±0.003% in less than one hour. Those 
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values separated by 0.050% which is about 12 times larger than the typical observational 
uncertainty in this sequence. There is also evidence of short-term variations in November 
5th, albeit with a smaller amplitude. For the July observations there is evidence of rapid 
polarization variations, but the case is not so clear because of the larger observational un- 
certainties. 

Time series analysis was performe d in the pola rization data using the Fast Fourier 
Transform with the Cleane st algo rithm (IFoster 1119951 ). Further details of the procedure can 
be found in IVinicius et al.l (120061 ). Two frequencies with high statistical significance were 
fo und: 0.074 and 0.57 c /d. The later frequency is close to the frequency of 0.49 c/d found 
by IVinicius et al.l (120061 ) . which they attribute to stellar rotation. This frequency, if real, is 
of great importance to our interpretation of the polarization variability in § |H 

In addition to variations in the polarization level, we also observed both short- and 
long-term variations in the PA. The long-term variations of the PA do not always correlate 
with the variations of the polarization level. For instance, from July 7th to 8th, we did 
not detect changes in the average polarization, but a 7° {6a) PA rotation was observed. In 
two-dimensional (2D), axisymmetric disks, the PA is expected to be constant and parallel 
to the symmetry axis. The observed changes are, thus, probably associated with departures 
from a simple 2D geometry. 



4. Modeling 



The above results raise the question of what are the physical mechanisms that cause the 
variations and whether the same mechanism is responsible for both the short- and long-term 
variations in the polarization. In the following we attempt to demonstrate that the observed 
polarization can be explained by the presence of blobs (small regions of density enhancement) 
and/or rings of gas close to the stellar photosphere. 

Before endeavoring to quantitatively explain the variations, we must find a suitable disk 
model that describes reasonably well the observed polarization and Ha line profile. For 



this purpose we use a recent three - dimen siona. 
described in ICarciofi fc BjorkmanI (120061 ) and 



radial disk density is con trolled by viscos ity (ILee et al 



non-local the r mody namic equilibrium code 



Carciofi et al.l (120061 ). We assume that the 



199ll ) and that the disk is hydrostatic 



in the vertical direction (jBiorkmanlll997l ). Given a prescription for the gas viscosity, the disk 
has two free parameters: the mass loss rate, M, and the outer radius, Rd- 

There is a wide scatter in the literature about the physical charac teristics of Achernar , 
since its rapid rotation makes it difficult to determine its spectral type (IVinicius et al.ll2006l ). 
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We adopt an spectral type of B3V, corresponding to a Tcfr ~ 20, 000 K and a radius of 
i?^ ~ 9Rq. We note, however, that varying those parameters within a reasonable range 
(say ±2000 K for T^s and ±li?0 for i?^) changes little the results below, a Eri has a 
strongly broadened Ha photospheric component and this must be acc ounted for when fittin g 



the observed profile. We have used the average spectrum of 1999 of IVinicius et al.l (120061 ) 



presumed to be purely photospheric, as an input photospheric profile in our modeling. 

In addition to the disk and stellar parameters, one must specify the angle i by which the 



system is observed. This parameter has been reasonably well constrained by IVinicius et al. 



(120061 ) to be around 70°. In our modeling we assume z as a free parameter in the range 60 



- 80°. 

We have chosen to model the observations of November 5th, since this is the only date for 
which we have simultaneous Ha and polarization observations. Our best-fit model is shown 
in Figure [3] and the adopted model parameters are summarized in Table [2l We fit exactly the 
5-band polarization level and our synthetic Ha profile reproduces well the observed profile. 

Having found a 2D model that reproduces well a subset of out observations, we now 
attempt to reproduce the polarization variations. Since this is an initial study, our goal is 
to find a qualitative explanation of the time dependence of the variations and a quantitative 
explanation of the amplitude of the variations. 

As stated in the introduction, the driving mechanisms operating in Be stars are still 
unknown. For this reason we adopt an ad hoc scenario which, nevertheless, has many of 
the physical ingredients believed to operate in Be stars. We assume that after a mass mass 
ejection event (hereafter MEE) a region of density enhancement, which we call a blob, is 
formed near the stellar photosphere. This blob is short-lived: after a few stellar periods (i.e., 
several days) the blob material will mix with the disk, both because of differential rotation 
and viscosity, thus forming a ring around the star. The fate of the ring is controlled by 
viscosity: if the mass transfer from the star stops (i.e., there are no subsequent MEEs), 
viscosity causes inner rim of disk to lose momentum and fall back to the star and the outer 
rim to diffuse outwards. If other MEEs occur, the ring density and size will increase with 
time. 



4.1. Rings 

Our assumption is that the observed long-term polarization variations are due to the 
formation and growth of a ring after one or more mass ejection events, or to the formation 
of an inner hole following a period of low mass loss rate. To model this, we define the 
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ring as a density enhancement in the disk between i?^ and Rr- The ring density is given by 
nr{r, z) = nd{r, z)'j, for r < Rj., where ra^ is the disk density and 7 is the density enhancement 
factor. An inner hole is represented 7 less than 1. 

In Figure m we show polarization results for models with Rr between 10 to 12 Rq and 7 
between and 3. From those models we learn that small rings with large density enhance- 
ments (7 > 2) successfully explain the observed amplitude of the long-term variations of the 
polarization. Similar effects are produced by larger rings with smaller density enhancements. 

We can estimate the time r required for the formation of a ring of radi us Rr using the 



following expression for the viscous diffusion timescale (e.g. iBjorkmaru 119971 ) 

T = Rliy-\ (1) 



where u is the kinematic viscosity of the gas. Assuming the eddy viscosity of lShakura fc Syunyaev 



(I1973I ) the above equation can be written as 

T = Rl {aaHy^ , (2) 

where a is the sound speed, a the viscosity parameter introduced by Shakura & Syunyaev 
and H is the disk scaleheight. For pressure supported, hydrostatic and isothermal disks, 
H can be written as H = aV~-^R~'^'^r^'^, where r is the distance to the star and Vent is 
the critical rotational velocity of the star. Assuming Vera ~ 400km s~^, a ~ 20km s~^ and 
a = 0.2, we find that a 2i?0-wide ring can be formed in a couple of weeks only. This means 
the timescale for ring formation or dissipation is of the same order of magnitude than our 
observed timescales for long-term polarization variations. 

We should point out that the ring model cannot account for the observed changes in 
the PA because it is intrinsically axisymmetric. In addition, it cannot explain the observed 
short-term variations, since the associated timescale of ring formation/dissipation (weeks) is 
much larger than the typical timescale of the short-term variations (hours) . 



4.2. Blobs 

We now attempt to quantify the polarization effects of an asymmetric gas distribution in 
the inner disk assuming a semi-spherical blob of material attached to the stellar photosphere 
at the equator. For simplicity we assume the blob has a constant density rj-biob- 

In Figure [5] we show the polarization level and PA as a function of the blob orbital 
phase for several blob sizes and densities. We adopt blob radii between 1 and 3Rq and ribiob 
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between 1 and 3 x 10^^ cm ^, which corresponds to 1.5 - 4.5 times the density at the base 
of the disk, hq. 

Our results show that blobs with the above characteristics can successfully explain the 
amplitude of observed short-term variations in the polarization level. For instance, the 
AP ^ 0.05% variation observed in September 14th can be accounted for by most models of 
Figure [H More importantly, however, the blob model can also reproduce the amplitude of 
the observed changes in PA. 

The blob model has two distinct associated timescales. The first is, clearly, the rota- 
tional period of the star, since in our assumption the blob is initially attache d to the stellar 



photo sphere. The orbital period of Achernar is very short (about 2 days, IVinicius et al. 



20061 ). and from Figure O we obtain that 6° changes in the PA or 0.05% changes in the 



polarization level can occur in A(f) < 0.1 or about 4 hours. 

The second timescale is associated with the process of the blob formation and growth 
during a MEE. It is difficult to estimate a timescale for this situation since we do not know 
how the blob forms; however, if we arbitrarily, for the sake of argument, give an initial 
velocity of a few hundred km/s to the blob material, we find that the blob can become as 
large as IRq in less than one hour. 

The above discussion is very qualitative and must be verified by detailed calcutations. 
The important point, however, is that the blob model appears to have associated timescales 
of one to a few hours, in accordance to our observations. 

The structures of the proposed scenario have a clear hierarchy: MEEs give rise to blobs 
which form rings which, in turn, form the disk. We cannot yet, with the present data, 
determine important quantities such as the characteristics (mass, geometry and rate) of the 
MEEs. However, from our modeling the typical masses of the structures involved can be 
roughly estimated. Our best-fit 2D disk model has a mass of 3 x 10~^° Mq. A ring of radius 
Rr = 11 Rq and 7 = 2, which can produce a polarization variation in the observed range, 
corresponds to a mass input to the disk of about 7% of the total disk mass. Similarly, a 
representative blob model, say with i?biob = 2 Rq and nbiob = 3 x 10^^ cm~^, corresponds to 
a mass input to the disk of about 5% of the disk mass. 



5. Conclusions 



We report the detection of both short- and long-term variations in the polarization level 
and position angle of Achernar. From an initial study using a state-of-the-art computer code 
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we arrived at a very plausible picture to explain the observed variations. In the proposed 
scenario, short-term variations are linked to two processes: the formation of a blob following a 
mass ejection event and the co-rotation of this blob with the (supposedly Keplerian) velocities 
of the inner disk. The finding of a significant frequency of 0.57 c/d in the polarization data, 
very close to the stellar rotation period, 0.49 c/d, strongly supports our assumption that the 
short-period variations are linked to photospheric processes. 

The long-term variations are probably a combinations of two processes: the formation of 
blobs and their subsequent dissipation to f orm a ring around the star. The ring scenario has 



already been proposed by iRivinius et al.l (120011 ) to explain line profiles variations observed 



in several early-type Be stars. Since this scenario is a natural consequence of the viscous 
decretion model combined with a non-constant mass loss rate, we suggest that observations 
aiming at detecting and characterizing rings in the inner disk of Be stars will constitute an 
important test of the viscous decretion model. 

The determination of the rate and typical masses of the mass ejection events would be 
a very important observational result. This cannot be achieved with our present data, but 
our results above strongly suggest that this determination should be possible with a long- 
term, high-precision polarization and spectroscopic monitoring of the early phase of the disk 
formation. 
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Table 1. Polarimetric Observations. 



Date Nr. of P PA 

(2006) WPPs (%) (°) 



Jul. 


07 


32 


0.145 ± 0.003 


35.8 ±0.6 


Jul. 


08 


24 


0.143 ±0.005 


29.1 ±1.1 


Jul. 


14 


8 


0.163 ±0.004 


31.2 ±0.4 


Jul. 


14 


20 


0.188 ±0.009 


29.5 ±1.3 


Jul. 


15 


16 


0.17±0.01 


34 ±2 


Sep. 


14 


34 


0.141 ±0.002 


30.6 ±0.4 


Nov. 


05 


41 


0.121 ±0.002 


34.6 ±0.5 
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Table 2. Adopted Stellar Parameters and Fitted Disk Parameters for Achernar 









M 


no 


Rd 


i 


(Re) 


(K) 


(kms~^) 


(Me yr-i) 


(cm~^) 




n 


9 


20, 000 


400 


8 X 10-^2 


6.7 X 10-^^ 


15 


65 
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Fig. 1.- 



B-filter polarization monitoring for a Eri. 
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Fig. 2. — Short-term polarization variability. Left column: results for September 14th. 
Right column: results for November 5th. The grey bands indicate the average value for the 
sequence, P, the band thickness spanning the value of la. 
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Fig. 3. 
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Best-fit model (thick lines) for the observations of November 5th. 
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Fig. 4. — Ring model. Shown is the B-band polarization vs. the ring outer radius, for several 
values of the density enhancement factor 7. The horizontal lines represent the minimum and 
maximum polarization observed by us. 
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Fig. 5. — Blob model. 5-band polarization (bottom) and position angle (top) vs. orbital 
phase. The models on the left show the effects of blob radius (-Rbiob = 1 to 3 i?©) for a given 
blob density (ribiob = 3 x lO^^cm"^). The models on the right show the effects of blob density 
('T'biob between 1 and 3 x lO^^cm"^) for a given blob radius (-Rbiob — 2.5 i?©). A phase of 
zero corresponds to the blob between the star and the observer. 



